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Abstract. The weak lensing effects are known to change 
only weakly the shape of the power spectrum of the Cos- 
mic Microwave Background (CMB) temperature fluctu- 
ations. I show here that they nonetheless induce specific 
non-Gaussian effects that can be detectable with the four- 
point correlation function of the CMB anisotropies. The 
magnitude and geometrical dependences of this correla- 
tion function are investigated in detail. It is thus found to 
scale as the square of the derivative of the two-point corre- 
lation function and as the angular correlation function of 
the gravitational displacement field. It also contains spe- 
cific dependences on the shape of the quadrangle formed 
by the four directions. 

When averaged at a given scale, the four-point func- 
tion, that identifies with the connected part of the fourth 
moment of the probability distribution function of the lo- 
cal filtered temperature, scales as the square of logarith- 
mic slope of its second moment, and as the variance of the 
gravitational magnification at the same angular scale. 

All these effects have been computed for specific cos- 
mological models. It is worth noting that, as the amplitude 
of the gravitational lens effects has a specific dependence 
on the cosmological parameters, the detection of the four- 
point correlation function could provide precious comple- 
mentary constraints to those brought by the temperature 
power spectrum. 

Key words: Cosmology: Dark Matter, Large-Scale Struc- 
tures, Gravitational Lensing, Cosmic Microwave Back- 
ground 



1. Introduction 

A robust prediction of inflationary scenarios is that the 
temperature fluctuations of the Cosmic Microwave Back- 
ground (CMB) are expected to obey Gaussian statistics. 
Actually this prediction has been challenged recently by 



several authors (Falk, Rangarajan & Srednicki 1993, Mun- 
shi, Souradeep & Starobinsky 1995) who calculated the 
skewness induced by nonlinear couplings in the primary 1 
stage of the temperature fluctuation generation. The skew- 
ness induced at this level has been found, however, to 
be entirely negligible compared to the cosmic variance, 
and thus not accessible to detections. Therefore, the pri- 
mary temperature maps are entirely defined, in a statisti- 
cal sense, by the power spectrum of the temperature fluc- 
tuations or equivalently by the shape of the two-point cor- 
relation function. A number of other statistical indicators 
are thus set up by this a priori hypothesis. In particu- 
lar Bond & Efstathiou (1987) have investigated expected 
properties of such temperature maps, as the number den- 
sity of temperature peaks, their correlation functions... 
Moreover, the exploration of the CMB physics has been 
boosted recently after it has been realized that it would 
be possible to determine all the cosmological parameters 
with a remarkable precision from an accurate, and acces- 
sible, measurement of the temperature power spectrum 
(Jungman et al. 1996). In particular, the effects of the 
secondary sources of temperature fluctuations (Sunyaev- 
Zcl'dovich effects, nonlinear Doppler effects, lenses..) and 
foregrounds (point sources, galactic dust) on the power 
spectrum have been investigated in more details (see Co- 
bras/Samba report, 1996, for a general discussion on these 
problems). These calculations have shown that most, if not 
all, of these effects have a relatively small impact on it. In 
these calculations, however, the impact of secondary ef- 
fects on the Gaussian nature of the temperature field has 
not been considered. Particularly interesting are the the 
higher order correlation functions that are identically zero 
for Gaussian fields, and are thus direct indicators of any, 
even small, non-Gaussian features. 

In this paper the calculations will be focused on the ef- 
fects of weak-lensing on CMB maps. They, indeed, consti- 
tute a particularly attractive mechanism because it comes 
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1 "primary" means the anisotropies induced on the last scat- 
tering surface by either potential fluctuations, Doppler effects 
or photon density fluctuations, not taking into account the sec- 
ondary effects nor the foregrounds. 
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from a coupling between the primary temperature fluctu- 
ation field and the mass concentration on the line of sight 
acting as deflectors. Their impact on the power spectrum 
has been investigated primarily by Blanchard & Schneider, 
(1987) who found the effect to be negligible. More recent 
works (Kashlinsky 1988, Cole & Efstathiou 1989, Sasaki 
1989, Tomita & Watanabe 1989, Linder 1990, Cayon, 
Martinez-Gonzalez & Sanz 1993a, b, Fukushigc, Makino 
& Ebisuzaki 1994, Seljak 1996) eventually confirmed this 
conclusion, although the point was debated for a while. 
In particular Seljak (1996) made a detailed calculations of 
these effects for realistic models of CMB anisotropies and 
using a semi-analytic calculation based on a power spec- 
trum approach that includes nonlinear corrections. In this 
text I will follow a rather similar approach to investigate 
the apparition of non-Gaussian features caused by weak 
lensing effects. 

In section 2, I present the basis of the physical mecha- 
nisms describing the CMB map deformations induced by 
weak lensing effects. In Section 3, I give the explicit ex- 
pression of the first non- vanishing correlation function, the 
four-point, in different remarkable geometries. In Section 
4, quantitative predictions are given for two different cos- 
mological models. The dependence of the results on the 
cosmological parameters, and the practical interests that 
such a measurement could have, are discussed in the last 
section. 



2. Weak Lensing Effects on CMB Maps 

2.1. The Basis of the Physical Mechanism 

The effect of a gravitational lens is to induce a displace- 
ment of the light path, thus moving the apparent position 
of a sky patch on the last scattering surface by a given 
angle. The temperature of this patch is not affected itself, 
i.e. lenses do not created new structures, and a perfectly 
isotropic sky would remain so. The patch of the sky ob- 
served at the position -y obs - is thus actually coming from 
the position -yP rim - on the "primordial sky" , and the dis- 
placement, 8-y, is induced by the mass concentration on 
the line of sight. More precisely 5~f is given by the trans- 
verse derivative of the projected potential <f> of the mass 
fluctuations, 
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where V is the angular distance, \ is the distance of the 
lenses along the line of sight and xcmb is the distance of 
the last scattering surface (see Kaiser 1992, Seljak 1996, 
Bcrnardeau, van Waerbeke & Mellier 1996 for more details 
on this equation). It is interesting to rewrite this equation 
in terms of the Fourier transform of the mass density fluc- 



tuation field. The Fourier transforms 6(k) are defined by, 

/c\ 3 k 
(2^)372 W X 

exp [iD fc± • 7 + ik r x] , (2) 

where the linear growth factor D + {\) and the Fourier 
transforms are normalized to the present time. Then the 
potential reads, 

,. , 3 D+(x) f d 3 fc 6(k) 

exp [iD k± ■ 7 + ik r x] (3) 
which implies that the displacement can be written, 
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with 



w(x) = 3^o 



A) (Xcmb - x) go(x) D +(x) 
£>o(xcmb) a(x) 



(5) 



The function w(x) gives the efficiency function of lenses 
for sources located on the last scattering surface. It will 
be investigated in more details in the last section. 

Note that in the following I will amply use the small 
angle approximation. It implies in particular that a given 
patch of the sky can be decomposed in flat waves and also 
that, in moment calculations, the component of k along 
the line of sight can be neglected compared to the norm 
of k±. 

2.2. The Effects on CMB Maps 

Compared to detections on background galaxies, the in- 
vestigation of lens effects on the last scattering surface 
is very attractive, because this surface is at a well de- 
fined rcdshift, and has a negligible width. The analysis of 
the lens effects requires however more sophisticated tools 
since the induced shear cannot be directly measured. The 
primordial temperature patches on the CMB sky arc in- 
deed known only statistically and have a large angular 
correlation length. In which way, then, can the lens effects 
be revealed? Actually lensed CMB maps can be seen as 
collections of temperature patches of different sizes and 
shapes, which or only a fraction of which are displaced 
or deformed. Although this is slightly arbitrary, two ef- 
fects can be distinguished in the way sizes and shapes of 
patches are affected, 

— the shear effect that deforms, stretches out tempera- 
ture patches in the shear direction, 

— the magnification effect that globally enlarges or 
shrinks those patches. 
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The local deformations of the temperature patches are 
however a priori difficult to disentangle from the actual 
primordial intrinsic temperature fluctuations 2 . What will 
make then the effects detectable is the fact that close 
patches will be deformed in a similar way (when they are 
seen through a unique lens), and the excess of these close 
rare features cannot be accounted from a Gaussian field. 
It is thus possible to quantify their presence by statistical 
indicators. The power spectrum is of course not adapted 
to take into account the apparition of such non-Gaussian 
features. For that matter the high-order correlation func- 
tions, that arc all identically zero for pure Gaussian fields, 
are extremely precious. Indeed these higher-order corre- 
lation functions contain informations about shapes, and 
their derivations can be pursued completely with Pertur- 
bation Theory techniques. In the following I focus my 
analysis on the first non vanishing correlation function, 
the four-point one. 

3. The Effects of Weak Lensing on Temperature 
Correlation Functions 

3.1. Statistical Properties 

From the first equation it is easy to see that the tempera- 
ture, <5y bs -, observed in the direction 7 is in fact the unaf- 
fected temperature #p" m ' coming from a slightly changed 
direction, 



S° bs -( 7 ohs -) = «5P rim -( 7 P rim -) = (5P rim -( 7 obs - + 5-y). 



(6) 



In the following I assume that the displacement is small 
compared to the angular scale at which the observations 
are made. This is a fair assumption since the displacement 
is at most of 1' (for cores of clusters) and that the angular 
resolution of the future satellite missions does not go below 
10'. As a result it is always possible to expand the relation 
(6) with respect to the displacement, 



^ bs '(7) 



^ lm '(7) 



d 
d~1 l 



(7)*7» 
(7)*7»*7j + • • • 



(7) 



where the Einstein index summation prescription is used. 

It is important to have in mind that both quantities 
$r" m (7) ana - ^7 are independent Gaussian fields. The pri- 
mordial temperature field is Gaussian for inflationary sce- 
nario, and independent of the lens potential field because 
the last scattering surface is far from the intervening lenses 
(that will be found to be at redshift ^ 5). The lens po- 
tential field will be assumed to be in the linear regime, 
although this is not a crucial hypothesis 3 . 

2 In particular there are no known working method to con- 
struct weak-lensing maps from CMB temperature maps. 

3 This assumption allows the use of the linear power spec- 
trum, but one could have included nonlinear corrections to 



In the following the correlation functions or moments 
will be calculated in the small angle approximation, for 
which the plane approximation for the last scattering sur- 
face can be made. Thus one can write, 

f d 2 l 

#r im '(7) = J j^) ai exp(i/- 7 ), (8) 

where the o; coefficients obey Gaussian statistical rules. 
In particular, 



(ai a t >) = <S D irac(* + C U 



(9) 



where the C; are the "famous" C; describing the angular 
power spectrum. 

On the other hand the random variables 8{k) obey the 
statistics, 



(S(k)6(k')) = 6 miac (k + k')P(k), 



(10) 



where P(k) is normalized to the present day. In order to 
produce a consistent set of power spectra it is important to 
have a consistent normalization for C; and P(k). This can 
be obtained from the small I behavior of C'i corresponding 
to the small k behavior of P(k). In the following I will 
assume an Harrison- Zel'dovich initial spectrum, so that, 

P(k)=Ak at small k. (11) 

The coefficient A can be related to the small I behavior of 
I (see Hu 1995 for instance), 

,(i + 1)Q = s n »(S)) 2 atmalu - ,12) 

3.2. The Two-Point Correlation Function 

The dominant corrective term for the angular two-point 
correlation function can be calculated from the expansion 
(7), 

« s -(7iX s -(72)) = (^ rim -(7i)^ rim -(72)) + (13) 
<^ rim -(7i)^(7i) ^ rim -(7 2 )^(7 2 )) + 

<^T im -(7i)^-^ rim -(7 2 )^(7 2 )^(7 2 )) + • • • 

The corrective terms have been written up to the 
quadratic term in the large-scale structure density field. 
The previous expression can be written in Fourier space, 

<<#»-(7i)$"-(72)> = Q exp[iZ. 7l2 ]+ (14) 

/ Cl exp [ iZ ^i2] J d Xidx 2 w(xi) w(x 2 )x 



P ( k ) ( ,2^ / ^ , ^ cxp[i(xi - X 2 )fc r ] x 



(2tt)3 v ' \k2V (xi)V (x2) 

(exp[iX> (xi)fc±-7i - iA)(x 2 )fc±'7 2 ] - 2) 



it, since the lens density field is by no means required to be 
Gaussian. 
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where 7 12 is the angular distance between 7l and 7 2 , 
7i2 = iTi — 1 ■ This formula has been obtained using eqs. 
(9,10) defining the two power spectra. To complete the 
calculations one can use the small-angle approximation, 
well verified below 1 degree scale, which implies that 



k^k ± and P(k)^P(k ± ). 



(15) 



Then the integral over k r leads to a Dirac function in 
Xi — X2- We eventually have, 



« s -(7iX s -(72)) = C( 7l2 ) + 



(16) 
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This result does not coincide apparently with the one of 
Seljak (1996, equation A. 6), but simply because the expo- 
nential was not expanded in his expression. That this ex- 
pansion can be done is amply justified by the fact that the 
displacements are small compared to the angular resolu- 
tion scale. Then one recovers exactly the same expression. 
Note that the autocorrelation function of the displacement 
field is automatically introduced by the third term in (13). 

3.3. Higher- Order Correlation Function 

It is quite easy to see that the weak-lensing effects do not 
introduce a three-point correlation function. It is indeed 
impossible to build a term of non-zero ensemble average 
involving three a/ factors. 

The first non trivial high order correlation function is 
thus the four-point correlation function. At this stage it 
is important to have in mind that the observable quan- 
tity is the connected part, 
(<5° bs - S° hs - (7 2 ) 6ft s - (73) 6ft s - (7 4 )) c , of the ensemble 

average, (<5° bs '(7i) ^ bs '(7 2 ) <^ bs (7 3 ) k^ili)), that is 
the part which is obtained when the products of two point 
correlation functions that can be built have been sub- 
tracted out, 



<^ bs -(7i; 



(7, 



(7a 



1 5° bs - 



(18) 



■(74)> c = 

(^ bs -( 7l )^ bs -(7 2 )^ bs -(7 3 )^ bs -(7 4 ))- 

<^ bs - (7i) ^ bs - (7 2 )) <^ bs - (73) <5t S - (74)) 
perm. (2 other terms). 



The connected part is obviously zero for the primordial 
field: it is a direct consequence of its Gaussian nature. 
The dominant term, in terms of weak lensing effects, is 
thus given by, 

( ( 5° bs -(7 1 )^ bs -(7 2 X s -(7 3 )'5T bs -(74)) c = (19) 



« s '(7i) ^-^ rim '(7 2 ))(^(7 2 ) 57,(73)) x 
(n^- <5 T im '(73) 5 T S '(li)) +pcrm. (11 other terms). 

7? 

Roughly speaking it means that the four-point correlation 
function, in units of the square of the second, is propor- 
tional to the weak lensing angular correlation function. 
Although at this stage it is difficult to give definitive quan- 
titative predictions, the magnitude of the fourth order cor- 
relation function should be about 10 -2 (the order of the 
corrective term in [13]), which should be easily detectable 
in full sky coverage CMB maps. The expression of the 
four-point correlation function can be given in terms of 
the power spectra, 



« s '(7i) ^ bs '(7 2 ) <St s -(7 3 ) ^t s -(7 4 )) c 



(20) 



d 2 h 



{2v) 2 C h J (2tt)2Q 
li ■ k li ■ k 
k 2 V k 2 V a 



dxw2{x) J Wf m 



exp [ili • 7 12 + iV k ■ 7 23 + iZ 2 • 7 34 ] 
perm. (11 other terms). 



This expression can be calculated with an integration over 
the angles between i x and 7 12 and l 2 and 7 34 respectively. 
It yields the Bessel functions </i(^i 7 i 2 ) and Ji(/ 2 7 34 ). 
The results can thus be expressed in terms of the angular 
derivative of the two-point correlation function, 



d 

d7 



C( 1 ) = -J l ^QJ 1 (l 1 ), (21) 

and with quantities associated with the angular correla- 
tion of the displacement field, 

D P (l) = J dx- 2 (x) J ^^§2 UVokj), (22) 
leading to (Appendix A), 

W ta -biH?''(7 2 )«?--(73)^ b ''(7 4 )) t = P3) 

[A)(7 23 ) cos(^i 2 - 1^34) - L> 2 (7 23 ) cos(i^i 2 + ip 34 )} + 
+ perm. (11 other terms), 

where ip\ 2 is the angle between 7 12 and 7 23 and (^34 is the 
angle between 7 43 and 7 32 (see Fig. 1). Two terms are thus 
involved. The a priori dominant term is the one in D , and 
it is weighted by the cosine of the angle ip\ 2 — (/? 34 = ip (see 
Fig. 1), that is the angle between the directions 7l2 and 
7 34 on the sky. It gives a clear geometrical dependence for 
the four point-correlation function. However, one should 
have in mind that 11 other terms have to be taken into 
account in this calculation. This signal may therefore be 
masked by other geometric dependences. 
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Fig. 1. Description of the angles intervening in the expression 
(24) of the four-point correlation function. The thick solid line 
materializes the dC(0)/d9 factor, whereas the hatched lines 
represents the correlation functions of the displacement field. 

Quantitative calculations can be done for specific cos- 
mological models (see next section). However, one obvious 
problem for a practical determination of this correlation 
function is that it depends on 5 different variables. It is 
thus crucial to reduce the number of variables by consid- 
ering simplified geometries. 



3.4- Peculiar Geometries 

3.4.1. When two Directions Coincide 

The first geometry one may think of is when two 
points arc merged together, that is the expression of 

(#r bs -(7i) [<$t s -(7 2 )] 2 ^t s, (7 3 )) c - This notation is actu- 
ally a bit oversimplified since the local temperature fluc- 
tuations are actually filtered by the used apparatus. One 
should thus have in mind that the two directions denoted 
7 2 arc actually close random directions in a beam centered 
on 7 2 . 

Of course, once again, many terms are contributing to 
this ensemble average but I will first concentrate on the 
case where the connection between the two 5^ >s ^-f 2 ) is 
made by the lens coupling term (subsection A.2.b). In such 
a case one can see that ip is given by the angle between 
72 — 7i an d 73 — 72 an d is not affected by the smoothing. 
This is not the case for the term in cos(<pi 2 + ^34) which 
is expected to vanish because it is averaged to zero (more 
precise derivations are given in the Appendix). This con- 
tribution is thus proportional to the cosine of the angle, 
and to the autocorrelation function of the displacement 
field. 

What about the other terms? Their geometrical repre- 
sentations are given in Fig. 2. The (2a) diagram is the term 
that has just been considered and is expected to dominate 
the final expression. Note that all these diagrams have a 
symmetry factor of 2 compared to what is given in the Ap- 
pendix. At first view the (2b) diagram vanishes because 
ipi2 takes a random value averaging both cos(<y9i2 — ^34) 
and cos((^i2 + 9334) to zero. A more detailed calculation is 
proposed in the Appendix. It shows that it gives a contri- 
bution proportional to the angular correlation function of 




Fig. 2. Geometrical representation of the terms intervening in 
the expression (25). 

the local weak lensing convergence (see Blandford et al. 
1991, Villumsen 1996, Bernardeau et al. 1996), 

Mo(7) = / dx^ 2 (x) J^P(k)MV k 7 ). (24) 

The two other diagrams are simpler since they can be 
simply obtained from the general expression (24). Taking 
all these terms into account we have, 

« s -( 7l ) [^ bs '(7 2 )] 2 <$T bs -(7 3 )) c « (25) 
^g^C (7l 3)[ 7l2 M 0(7l2 )cos(^ l) + 

723 M o(723) C0S(V> 3 )] + 

^C(7i 2 ) ^ C '(723) C0S W x 

[D (6 a ) + A}(7i 3 ) - A>(7i 2 ) - A)(7 23 )] 

where ip 3 is shown on Fig. (2b) (tpi corresponds to the di- 
agram obtained when the roles of 7l and 7s are inverted) 
and 9q is the smoothing angle of the experiment. In this 
expression, the D2 terms have been neglected, Do has been 
changed in Do for diagram (2a) and C in C in (2b). This 
is due to the filtering effects. One should indeed take into 
account the average displacement within the beam size 
of the experiment in the first case, the rms temperature 
fluctuations in the second. More precisely we have, 

Do(0) = J d xw *( X ) J*£ |^ W\V o k0), (26) 

and 

C(0)= J l ^C,W 2 (l6), (27) 
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with, for a top-hat window function for instance, 

W {x) = (28) 



It is then interesting to define the function, 

«4 (712,723)= ( 29 ) 
(g^Ti) [^ bs -(7 2 )] 2 ^T bs -(7 3 )) c 
(^ bs -(7i)^ bs '(7 2 )) <^ bs -(7 2 X s -(7 3 )) ' 

which is a dimcnsionlcss quantity. It does not depend in 
particular on the magnitude of the CMB temperature fluc- 
tuations and it is directly proportional to the large-scale 
structure power spectrum, with a known dependence on 
the shape of the anisotropy power spectrum. This quan- 
tity would thus be a measure the weak lensing effects in 
CMB maps. Taking into account the fact that the diagram 
(2b) generally leads to a negligible contribution we have, 

«4(7l2>723)~ (3°) 

^log[C(7i 2 )] ^log[C(7 23 )] cos(^) x 

[Aj(0o) + £0(713) - £0(712) - £0(723)] • 

A 2D contour plot of the function «4 (712)723) is proposed 
in Fig. 6 for a peculiar cosmological case. 

The main contribution to this expression is the term 
coming from the diagram (2a), but the contributions from 
(2c) and (2d) cannot be neglected because the correlation 
function of the displacement field £0 is only slowly de- 
creasing with the angle (see next section) . It is interesting 
to have in mind the physical effect described by this re- 
sult. It corresponds indeed to a shear effect that is the de- 
formation of two nearby temperature patches by a single 
lens. It can thus be seen as the excess of close temper- 
ature peaks that are elongated in the same direction. It 
is therefore logical that this effect is proportional to the 
correlation function of the displacement field. However, 
it does not always dominated the four point correlation 
function. This is the case in particular at the degree scale 
since the derivative of the temperature correlation func- 
tion vanishes (see Fig. 4 in the next Section). The next 
subsection is devoted to a simpler geometry where this 
case is more specifically investigated. 

3.4.2. When two Pairs Coincide 

This case is obtained from the previous case when the 
directions 7 3 and 74 coincide. The geometrical represen- 
tations of the involved terms are presented in Fig. 3. The 
a priori dominant term corresponds to the diagram (3b). 
However, as noted previously, both the diagrams (3b) and 
(3c) vanish for specific values of the distance between ■y 1 
and 7 2 . In such a case we are left with the diagram (3a). 
A crude evaluation of this diagram gives zero because the 



1 Y " (3a) 




(3b) 



(3c) 



Fig. 3. Geometrical representation of the terms intervening in 
the expression (31). 



averages over the angles <^i2 and (^34 are expected to van- 
ish. However this is true only when the smoothing angle 
is negligible compared to the angular distance between j x 
and 7 2 . More precise calculations (subsection A. 3) show 
that this term is proportional to the logarithmic derivative 
of the temperature variance in the beam. The resulting 
value of the correlation function is (taking into account a 
symmetry factor of 4), 



U^ bs '(7i)] [^ bs -(7 2 )] > 
2 



(31) 



1 2 



[£ (^0) - £0(712)] • 

Note that the first term is no more proportional to the 
angular correlation of the displacement field, but to the 
angular correlation of the convergence (that is, within a 
factor 2, of the magnification in the weak lensing regime). 
Here the physical mechanism has changed. The correla- 
tion function is not due to a local shear, but to single 
lens amplification of patches on the primordial sky that 
create an excess of close bright peaks. It is thus clearly a 
magnification effect. 

3.4.3. When the four Directions Coincide 

The last case I consider is when the four directions coin- 
cide. In this case the magnification effect always dominates 
and the results, obtained in subsection A. 4, are similar 
to the one discussed in the previous subsection, but with 
small changes introduced by the filtering effects, 

(P'W] 4 ),^^) 2 ^^) (32) 



with 



M O (0 O )= / d X w 2 ( X ) 



fcdfc 
~2tT 



P(k)W 2 (V o k0 o ). (33) 
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This result has been properly demonstrated for a top-hat 
window function, and should be roughly correct for other 
window functions. Note that this expression is also the 
fourth cumulant of the local filtered temperature prob- 
ability distribution function. We can as well define the 
dimensionless quantity, 



([S^-iOo)] 4 



>c 



■(°o)] 2 ) 2 



dlQg[C(M \ ,« M n A , 

dlog[0] ) Mo ^ 3 ^ 



which tells us that the dimensionless kurtosis of the local 
filtered CMB temperature probability distribution func- 
tion is proportional to the variance of the local filtered 
convergence. 

4. Quantitative Predictions 

4-1. The Cosmological Models 

The cosmological models used to illustrate the previous re- 
sults by quantitative predictions are standard CDM mod- 
els with fibaryon = 0.05, Hq — 50 fcm/s/Mpc with an ini- 
tial Harrison-Zel'dovich spectrum. Two cases have been 
chosen, f2 = 1, A = (model 1) and Q = 0.3, A = 0.7 
(model 2). The transfer function and temperature power 
spectrum were both computed with the code of Seljak & 
Zaldarriaga (1996). 

For convenience the mass fluctuation power spectra 
were approximated by simple analytic fits (similar to the 
ones proposed by Bond & Efstathiou 1984), 



provide the normalization constraint for A, 
A w 1.01 HT 8 O 2 ( Q( ? CMB) - V. 

\D+(XCMB)J 



(39) 



Note that in the previous equations, distance units have 
been chosen so that cH = 1. 

4-2. Numerical Results 



CD 

u 



100 




0/deg 

Fig. 4. The function 1/C(0) dC(9)/d8 as a function of the 
angle. The thick line is for model 1 and the thin line for model 
2 



P(k) = A 



(1+ [afc+(6fc) 3 /2 + (cfc) 2 ] tl ) 2/l ' 



with 



a = 



1.13; 
6.5 

3000 r ' 

3 

3000 r ' 

1.7 

3000 r' 



and 



r = 0.5 for ft = 1.0 
r = 0.13 for Q a = 0.3. 

The normalization factor is given by the relation (12) with 
observational constraints, 



1(1 + I) a = 



24tt fQo 



The measured values of Qo and T by the COBE satellite 
(Mather et al. 1994, Gorski et al. 1994) 

T = 2.726 ± 0.010 K, 
Qo = 19.9 ± 1.6 10~ 6 K, 



(35) I present the derivatives of the temperature angular 
correlation functions in Fig. 4. They exhibit a remarkable 
property, since they both drop to zero at a scale slightly 
larger than 1 degree. This property can be of great help 
to disentangle various contributions. For instance it im- 
plies that specific geometries produce a vanishing four- 
point correlation function. This is the case for instance 
when three of the four directions form an equilateral tri- 
angle with ~ 1.2 deg of side length. Such a property could 
be of crucial interest to ascertain the origin of an observed 
four-point correlation function. 

The two other quantities of interest are the functions 
A) (7), D 2 {i) and A) (7) (eqs. 22, 26) that describe the 
magnitude of the lens effects. They are presented in Fig. 
5. One can see that Do and Dq dominate at small scale. 
The resulting shape of the dimensionless four-point corre- 
lation function «4(7 12 , 723) (eq. 29) is presented in Fig. 6. 
It exhibits specific features induced by the cos(V>) factor 
and by the derivative of the temperature angular correla- 

(36) tion function. The fact that the latter vanishes is clearly 
present with a significant circular feature at 1 degree scale. 
It implies that the kurtosis is maximum for angular dis- 
tances below 1 degree. Note that in this figure, all con- 
tributions have been included, but the contributions from 

(37) (2b) are found to be negligible and not to affect the global 

(38) features of this plot. 
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0/deg 

Fig. 5. The functions Do (solid lines), D2 (dashed lines) and 
-Do (long dashed lines) for model 1 (thick lines) and model 2 
(thin lines). 
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U 0.4 
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6/deg 

Fig. 7. The functions dlog[C(0)]/dlog[0] for model 1 
line) and model 2 (thin line) 
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Fig. 6. Contour plot of the function K4(7 12 , 723) ( ec L- 29) as 
a function of the relative position (in degrees) of 7 3 when 7 2 
is the central point of the graph and 7 X is at the coordinates 
x = 0, y = +0.63 deg. The value of K4 has been multiplied by 
1000. 



I also present the quantities intervening in the expres- 
sion of the dimensionless kurtosis of the CMB temperature 
PDF. The index of the local temperature fluctuations is 
given in Fig. 7 for the two cosmological models. The vari- 
ance of the convergence is presented in Fig. 8, and the 
resulting value of K4 in Fig. 9. 

The resulting coefficient K4 can be as high as 5 10~ 3 . 
It depends however a lot on the cosmological models 
and more particularly on the amount of power at about 
10ft. _1 Mpc scale. The reason why Mq, and consequently 
K4, is smaller for the fl = 0.3 case is thus actually due to 
the change of shape of P(k) (i.e., a lower value of T), and 
not directly to the low value of f2. The dependence on the 
cosmological parameters is investigated in more details in 
the next subsection. 



IS 

a 




0/deg 

Fig. 8. The functions Mo (solid lines) and Mo (dashed lines) 
for model 1 (thick lines) and model 2 (thin lines) 




0.1 



Fig. 9. The function ka{0) 
and model 2 (thin line) 



0/deg 

(eq. 34) for model 1 (thick line) 
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4-3. Dependence on the Cosmological Parameters 

The dependence on the cosmological parameters enters in 
various aspects. It is important in particular for the rela- 
tion between the C; and the mass fluctuation power spec- 
trum P(k). This is not the aim of this paper to explore 
in great details all these dependences. There is however 
a dependence which is specific of the lens effects, that is 
the shape of the efficiency function w{z) (cq. 5). In Fig. 
10, I present the function w(z) for different cosmological 
models. We can see that both the shape and the ampli- 
tude of this function strongly depend on the cosmological 
parameters. 

The main dependence is actually due to the factor Q, 
appearing in the expression of w(z). This is however quite 
misleading since this factor cancels out with the one in- 
tervening in the expression of the normalization factor A 
(eq. 39). The dependence on the cosmological parameters 
is thus mainly contained in the shape of the power spec- 
trum, that is in the ratio between the amount of power 
at relatively small scale to the power at very large scale, 
although they enter significantly in the efficiency function 
w(z), even when the overall fl factor has been dropped. 
The way the two cosmological parameters f2 and A enter 
in this expression is unfortunately rather cumbersome. It 
is however clear that the CMB four-point correlation func- 
tion contains informations on the cosmological parameters 
in a quite different combination compared to the one inter- 
vening in the temperature power spectrum. In particular 
the detection of such an affect may allow to better disen- 
tangle what is specific of C; from what contribute directly 
to P(jfe). 




10 

redshift z 

Fig. 10. The function w(z) (eq. 5) for the Einstein-de Sitter 
case (thick solid line), a model with £7 = 0.3, A = 0.7 (thin 
solid line) and a model with Q = 0.3, A = (dashed line). 



5. Discussion and Conclusions 

In this work, I have calculated the expression of the 
four-point temperature correlation function as induced by 



weak-lensing effects. For standard CDM model, the ampli- 
tude of this correlation function, in units of the square of 
the second, is found to be of order of 5 10~ 3 . However, this 
estimation did not take into account the nonlinear evolu- 
tion of the power spectrum that might significantly am- 
plify this signal at small angular scale. This is for instance 
what is predicted for the two-point correlation function 
of the polarization of background galaxies (Jain & Sel- 
jak 1996). Unfortunately, in the case of CMB maps, the 
scale at which this effect might appear cannot be deduced 
straightforwardly from this work. This effect is indeed the 
result of a line of sight integration that mix different scales 
and different redshifts for a given selection function (e.g. 
Fig. 10) which is itself dependent on the redshift of the 
sources. Moreover all the intervening quantities have non 
trivial dependences on the cosmological parameters that 
should be taken into account. A detailed examination of 
the nonlinear effects is then left for a forthcoming paper. 

I would like to stress, that the amplitude of the lensing 
effects should be large enough to be detectable, at least 
marginally, in full sky CMB anisotropy measurements. 
The possibility of doing such measurements is directly re- 
lated to the cosmic noise associated with the quantities 
of interest, fourth moment or four-point correlation func- 
tion. So far, no precise estimation of the cosmic noise for 
the four-point correlation function has been made, but 
following Srcdnicki (1993), who presented the calculation 
for the three-point correlation function, one expects the 
cosmic variance of those quantities to be of the order of 
1/^pk. where l p k. is the typical value of I contributing to 
the temperature fluctuations. One can observe that the 
value of 1/Zpk. is of the order of the signal, however, one 
should have in mind that a direct and too naive calcula- 
tion of the cosmic noise may be actually misleading since 
the long wavelength fluctuations (corresponding to the low 
I part of the power spectrum) contribute significantly to 
it, whereas the lensing signal originates mainly from the 
small angular scales (below 1 degree). It suggests that 
the weak lensing signal might be more easily detectable 
in maps where the long wavelength temperature fluctua- 
tions have been removed. Moreover the detailed methods 
used to extract the signal might also be of different ro- 
bustness against the cosmic noise. In particular, it could 
be fruitful to take advantage of the a priori knowledge of 
the geometrical dependence of the four-point correlation 
function (see for example the cos(ip) factor in the expres- 
sion [30] ) . In a forthcoming paper, we explore the different 
possible strategies for the data analysis, and will present 
detailed estimations of the precision at which such a de- 
tection could be made in the future satellite missions. 

It also has to be noted that other secondary effects or 
foregrounds may also contribute to the four-point corre- 
lation function, not to mention the case of more exotic 
cosmological models based on intrinsically non-Gaussian 
topological defects. In particular the nonlinear Dopplcr ef- 
fects could induce a significant four-point correlation func- 
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tion, because it is caused by intrinsically non-Gaussian 
objects. There are however no reasons for these effects to 
have the same geometrical dependences, cos(ip) factor and 
dependence on the shape of the temperature two-point 
correlation function. Hence, it should be possible to dis- 
tinguish this effect from other sources. 

The most exciting aspect of this analysis is probably 
that the magnitude of the effect depends on the cosmo- 
logical parameters, 0, A and P(k) in a known way. The 
detection of the temperature four-point correlation func- 
tion may thus reveal to be extremely precious to test the 
global picture of the large-scale structure formation, as it 
will be unveiled by CMB anisotropy measurements. 
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Appendix A: Calculation of the Filtering Effects 

In this appendix I explicitly take into account the filtering 
effects to compute the expressions of the four-point cor- 
relation function in different geometries. The filtering can 
be the one due to the angular resolution of the apparatus 
or due to a subsequent filtering of the temperature field. 
When it needs to be specified the adopted window func- 
tion will always be the angular top-hat window function. 

The expression I am interested in is the expression (19) 
where -yP rim - and the coupling term V7 prim ' • £7 are given 

by 



W (7)= / ^a lW (l9 )e^, 
and, 



(40) 



V7 prim -(7) • £7 = (41) 

/(0^w/dxMx)/70x 

ai S W T^T W (\ l + V ° k \ o) e MKl + A)*0 • 7], 

where is the smoothing angle. The quantity of interest 
is thus 

<?4(7i,7 2 ,73,7 4 ) = / ^2 C h / ^f Cl * X ( 42 ) 



h k 

k 2 V Q 



<l\«- 2 (\) / ^P(k)W(h6o)x 



exp[iii • 7 12 ] W (\h - T> k\ ) exp[iD fc • 723] 

U k 



W (|/ 2 - V k\ Oo) exp[ii 2 • 7 34 ] W(l 2 9 ). 

In the following I will estimate this expression for different 
hypothesis on ■y 1 , 7 2 , 7 3 and 7 4 . 



A. 1 Four separate Directions 

Here I assume that 

0o <C 7jj whatever i and j. 



(43) 



The integral (42) will be dominated by values of l±, I2 and 
k for which /i7 12 , ^2734 and kj 23 are about unity. It 
implies that k 9o and |f , — T>ok\ 9o are all small quantities 
thus making the filtering effects negligible, so that, 

C4(7l,72,73,74)~/ (^2^ / ^f Cl * X ( 44 ) 

J d X w 2 ( X ) j J^ p (k)x 
exp[iZi • 7 12 + \V k ■ 7 23 + il 2 ■ 7 34 J 



k 2 V 



k 2 Vo 
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2 1. 

x 



An interesting property to be used is that J d\w 2 (x) J - - P(k) 

/ (^ -^-^It^ 1 ^ exp[iP 0fc • 7l2 + i* 2 • 723 ] ^ x 

It implies that (V 2 ^ O ) ^ - W'(h 9 ) W(h 9 ) { \^P) • 

C 4 (7i, 7 2 , 7a,7 4 )= (46) 

dC(7) dC( 7 ) f When one integrates over the angle of l\ the first term 

~~ d~7 — ^ 12 ^ — d-y — ^ 34 ) J dxw 2 (x) x vanishes. The second term of the expansion is thus the 



dominant contribution, which takes the form, 



f tfk_ _P(fc[ 712'fc 734 •*! ,~ , 

J {2Kfk*Vl 7l2 fc 734 fc ex P^ofc-7 23J , ^ ^ ^ L^irW,) / >R, 

C(7)= IS^Ci^. (47) 



with, 



^4(71,71,72,73) = -4^0— (2^)2 c ' 2 (55) 

I d X ™ 2 (x) I 7^ P ( fc ) x 



exp[iX> fc • 7i2 + ll 2 ■ 723] 

This integral can be eventually integrated from properties 

of Bessel functions, Using the property (45) we have 

/ d0exp[ia;cos(0)] cos 2 (^) = tt [J (ar) - J 2 (a:)] , (48) C 4 (7i,7i,72,7 3 ) = -^0-^- ^C( 723 ) x (56) 
/ d0 exp[ia;cos(t/>)] sin 2 (</>) = tt [J (x) + J 2 (x)] , (49) 7l2 ' 723 / d X w 2 ( X ) x 

Jo 723 J 

fcdfc J (P ^7i2) + ^2(^0^712) 



/■27T f kdk 

J d<j> exp[ixcos((/>)] sin(0) cos(^) = 0, (50) y -7^7 p ( fc ) 



with which we find, Interestingly C4 is now proportional to the angular corre- 

lation function of the local magnification and not of the 
(6° bs -( 7l ) <S° bs -( 72 ) <^ bs '(7 3 ) <$T bs, (7 4 )) c = (51) local displacement. 

2~frf C{ni2) fr^ C{risi)X A.2.b When 7 2=73 

[A>(7 23 ) cos((^ 12 - <^ 34 ) - £2(723) cos(<^ 12 + ^34)], In this case k is expected to be of the order of 1/0 O , thus 

w ^ ere larger than ?i and Z 2 . As a result one has, 

A>(7) = /dx- 2 (x) f W^T), (52) ^(7i,72,72,7 3 ) = / 7^2^/(1^ ^ (57) 

/2 /" d 2 fc 
dx ™ (x) y (2^ p(fc) x 

between 743 and 7 32 (see Fig. 1). f fe f fe 

exp[iZi • 7l2 + iZ 2 • 7 23 ] tJ^ - w2 ( k 9 ° ) tJ^ - ' 

.4. ,2 When two Directions Coincide 

A.2.aWhen 7l =7 2 leading to 

In this case, Zi is expected to be of the order of 1/6 Q , thus C 4 (7 1 ,7 2 ,7 2 ,7 3 ) = -— C*(7 12 ) — C(7 23 ) x (58) 
larger than k so that 2 d7 d7 

cos(V>)A)(0o), 

W(\h-V k\9 )^W(he )- (53) 

_^*fl 2 ?bW '(I 1 £b) + ... Wh6re 

As a result one has, ^ = / 1ST ^ < M > 

f d 2 li f d 2 Z 2 Here the effect is proportional to the mean displacement 

M7i,7i,7 2 ,7 3 )=J^(' il j— G l2 (54) in thc beam size ^ 
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A.2.C Other Cases To complete the calculation it is interesting to have in 

mind the property (Bernardeau 1995), 

The other cases do not give specific formulae and can be 
derived from results of section A.l. 



/lw^+<4 + ¥)= (64) 



A. 3 When two Pairs Coincide f hdh Tlr/ , . f hdh f_„ n 1 



W{h) + -h W'(l 2 ) 



2tt 

A. 3. a When 7j = 7 2 and 7 3 = 7 4 

_,, . .... , . _ . , , , . This property is rigorously exact. In principle it is not pos- 

Ihis case is similar to the A.z.o case where the result is ., , , , ,, „ .. , ,. n 

. . . . . . . . . . sible to separate the two terms. Both relations are however 

dominated by the second order term ol an expansion m , . ., .... . . , . . 

, ... . ,. , , ,,. ,. good approximation as it will be shown m the following. 

k0 , that here should be written for l\ and t 2 . We then T , . , , 

. It is thus reasonable to assume that 

have 

*7,,^-/^/^<i fg?w w * + " h r~ (65) 

f^tof^m* (») y^iw"^ 

//, . k\ 2 Then, using this property it is easy to show that 
l 1 9 W(l 1 9o)W'(l 1 9 ) (jj-j exp[iP fe • 7i 2 ] x 

^)\ w ^<<w *-_( s S9P) , **»- w > 

The integrations over the angles between h and k and A. 5 Validity of the top-hat window function property (65) 

between Z 2 and k give each a factor 1/2 leading to _ ,. r > r 

lo examine tnc property (bo), an interesting property of 

_ 2 the Bessel function to use is that (Gradshteyn & Ryzhik, 

C4(7i,7i,7 a ,7 9 ) = Y6 X < 61 > 198 °' ^ [8 ' 532 - 1]) ' 



[a 2 ( , f d 2 k 2J 1 (|fc 1 +fc 2 |) — v 2Ji(fei) 2J 1 (fc 2 ) 

J dxW (X) 7 J^P(k)eM&ok. ll2 ], ______ __ 

which can be expressed in terms of the angular correlation C 1 f — ^ 2 ] 

function of the magnification. V k\kz J 



(67) 



with 



A.3.b When •y 1 = 7 4 and 7 2 = 7 3 

This case is a pa rti c„,„ case of subsection A.2.a. C l I —<«"]= W <«) 



We have thus 

f-271 



j4.^ WTien i/ie four Directions Coincide 

In this case we have / dip Cp(~ cos(ip)) — for p odd, 

Jo 



(69) 



/d k lo ■ k U ■ k f 

(2^ PW F%F% X (62) I d^(-cos(<,)) = _r for p even. (70) 

W (|/i - P fe| O ) W {\h - V k\ O ) W(h 9 ) W(h 9 ). 

As a result 

This expression cannot be simplified furthermore if the /■ ^ /" dfc /" dZ 2 

window function is not specified. In this paragraph I as- ~ 2 — Ci 1 j — P(k) / — - C/ 2 (71) 

sume that W is the top-hat window function, _. 

Ji(he ) Ji{he ) (i + 2pi)(i + 2p2) x 

2 Ji(x) PuP2 
W{x) = ——. (63) J 1+2pi (/^ ) J 1+2pi (fc0 o ) Ji+2 P2 (fce ) Ji+ 2 p 2 (/ 2 0o) 
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In the following I assume a power law behavior for both 
the mass power spectrum and Cz, 



P(k) oc k nm , 
Ci oc 

Then using the property, 



r(A) 

2 X 



J„(t)Mt)t- x dt 



p | v+fj, — X+l 



(72) 
(73) 



(74) 



t/+M+A+l 



-v+fi+X+1 



one can easily compute the few terms of the previous se- 
ries. It is found to converge very rapidly. The resulting 
ratio, 



r(n c , n w ) = J^ Pl , P2 (1 + 2pi) (1 + 2p 2 ) x 
S(pi, l,n c ) S(p 1 ,p 2 ,n w ) S(l,p 2 ,n c ) 



(75) 



5(l,l,n c ) 2 5(l, l,n w ) 



with 



(76) 



S{Pi,P2,n) =— — -— X 

r(p 2 -pi+i-«/ 2 ) 
i 

T (p 2 + Pl + 2-n/2) r(-p2+ Pl + l-n/2)' 



is plotting in Fig. 11. It shows that the error made by 
using the approximation is at most of a few percent for 
the values of n r and n,„ of interest. 
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Fig. 11. The ratio (75) as a function of n c and n u 
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